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We present evidence for D°-D° mixing in D° K+Ti- decays from 384 fb"^ of e+e" coUiding- 
beam data recorded near = 10.6 GeV with the BABAR detector at the PEP-II storage rings 
at SLAG. We find the mixing parameters x'^ = [-0.22 ± 0.30 (stat.) ± 0.21 (syst.)] x 10"^ and 
y' = [9.7 ± 4.4 (stat.) ± 3.1 (syst.)] x 10~^, and a correlation between them of —0.94. This result 
is inconsistent with the no-mixing hypothesis with a significance of 3.9 standard deviations. We 
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measure Ru, the ratio of doubly Cabibbo-suppressed to Cabibbo-favored decay rates, to be [0.303 ± 
0.016 (stat.) ± 0.010 (syst.)]%. We find no evidence for CP violation. 



PACS numbers: 13.25.Ft, 12.15.Ff, ll.30.Er 



We present evidence for charm-meson {D^-D'^) mix- 
ing. This work complements results in the neutral 
-ft^ [H, H, -B i, d, and J5s [1, i] systems. Although 
precise predictions are difficult, D^-D^ mixing in the 
Standard Model (SM) is expected at the 1% level or 
less 0, B S ES El, [i3- Our result is consistent 
with this expectation and previous experimental lim- 



181. By observing the wrong-sign 
19|, we determine -Rd, the ratio 



its [H 0, 0, H [T? 

decay D° K+tt' 
of doubly Cabibbo-suppressed to Cabibbo-favored decay 
rates, and the mixing parameters x'^ and y' . 

The charm sector is the only place where the contribu- 
tions to CP violation of down-type quarks in the mixing 
diagram can be explored. We compare results for D'~' 
and decays but find no evidence for CP violation. 
The SM predicts effects well below the sensitivity of this 
experiment. 

We study the right-sign (RS), Cabibbo-favored (CF) 
decay — + K~7t~^ and the wrong-sign (WS) decay 
Z?" K^TT^ . The latter can be produced via the doubly 
Cabibbo-suppressed (DCS) decay K^tt^ or via 

mixing followed by a CF decay ^ ^ K^n^ . The 
DCS decay has a small rate i?D of order tan"* 9c ~ 0.3% 
relative to CF decay. We distinguish and by their 
production in the decay _D*+ "'^-D'' where the 7r+ 
is referred to as the "slow pion" . In RS decays the 7r+ 
and kaon have opposite charges, while in WS decays the 
charges are the same. The time dependence of the WS 
decay rate is used to separate the contributions of DCS 
decays from D^-D^ mixing. 

The and D° mesons are produced as flavor eigen- 
states, but evolve and decay as mixtures of the eigen- 
states Di and D2 of the Hamiltonian, with masses and 
widths Afi, Fi and M2, F2, respectively. Mixing is char- 
acterized by the mass and lifetime differences AAf = 
Ml — M2 and AF = Fi — F2. Defining the parameters 
X = AAf/F and y = AF/2F, where F = (Fi -I- F2)/2, we 
approximate the time dependence of the WS decay of a 
meson produced as a D'^ at time t = in the limit of 
small mixing (|a;|, \y\ <C 1) and CP conservation as 



TwsW 



-Tt 



oc i?D + Vi?D2/' Ft 



y 



.(Ft)2, (1) 



where x' = xcos5k-k + ysmSK-n, y' = —xsmSK-n + 
y cos 5 K-K, and 5k-k is the strong phase between the DCS 
and CF amplitudes. 

We study both CP-conserving and CP-violating cases. 
For the CP-conserving case, we fit for the parameters 
Pd, a;'^, and y' . To search for CP violation, we apply 
Eq. ^ to and samples separately, fitting for the 



parameters {Pq, x''^^ , y'^} for (+) decays and D'^ 
{—) decays. 

We use 384 fb~^ of e+e" colliding-beam data recorded 
near ^/s = 10.6 GeV with the BaBAR detector ^ at 
the PEP-II asymmetric-energy storage rings. We se- 
lect -D° candidates by pairing oppositely-charged tracks 
with a K^n^ invariant mass niKTr between 1.81 and 
1.92 GeV/c^, requiring each track to have at least 12 co- 
ordinates in the drift chamber (DCH). Each pair is iden- 
tified as K^n^ using a likelihood-based particle identifi- 
cation algorithm. The identification efficiency for kaons 
(pions) is about 85% (95%); the misidentification rate of 
kaons (pions) as pions (kaons) is about 2% (6%). 

To obtain the proper decay time t and its error 6t 
for each candidate, we refit the and tt^ tracks, 
constraining them to originate from a common vertex. 
We also require the D'^ and 7r+ to originate from a 
common vertex, constrained by the position and size 
of the e^e~ interaction region. We require the 7r+ to 
have a momentum in the laboratory frame greater than 
O.lGeV/c and in the e+e" center-of-mass (CM) frame 
below 0.45GeV/c. We require the probability of the 
vertex-constrained combined fit P(x^) to be at least 
0.1%, and the mD,+ —mKir mass difference Am to satisfy 
0.14 < Am < 0.16GeV/c2. 

To remove candidates from P-meson decays and 
to reduce combinatorial backgrounds, we require each 
to have a momentum in the CM frame greater than 
2.5 GeV/c. We require -2 < i < 4ps and 5t < 0.5 ps (the 
most probable value of St for signal events is 0.16 ps). For 
D*"*" candidates sharing one or more tracks with other 
D*~^ candidates, we retain only the candidate with the 
highest P(x^). After applying all criteria, we keep ap- 
proximately 1,229,000 RS and 64,000 WS D" and D" 
candidates. To avoid potential bias, we finalized our data 
selection criteria and the procedures for fitting and ex- 
tracting the statistical limits without examining the mix- 
ing results. 

The mixing parameters are determined in an unbinned, 
extended maximum-likelihood fit to the RS and WS data 
samples over the four observables mK-m Am, i, and 5t. 
The fit is performed in several stages. First, RS and WS 
signal and background shape parameters are determined 
from a fit to mK-n and Am, and are not varied in subse- 
quent fits. Next, the proper-time resolution function 
and lifetime are determined in a fit to the RS data us- 
ing rriK-n and Am to separate the signal and background 
components. We fit to the WS data sample using three 
different models. The first model assumes both CP con- 
servation and the absence of mixing, and only measures 
Pd. The second model allows for mixing, but assumes no 
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CP violation, and the third model allows for both mixing 
and CP violation. 

The RS and WS {rriKTr, Am} distributions are de- 
scribed by four components: signal, random 7r+, mis- 
reconstructed D° and combinatorial background. Signal 
has a characteristic peak in both mK-w and Am. The 
random 7r+ component models reconstructed decays 
combined with a random slow pion and has the same 
shape in mKjr as signal events, but docs not peak in 
Am. Misreconstructed events have one or more of 
the D'^ decay products either not reconstructed or re- 
constructed with the wrong particle hypothesis. They 
peak in Am, but not in rriK-K- For RS events, most of 
these are semileptonic decays K~t^v with the 

charged lepton misidentified as a pion. For WS events, 
the main contributor is RS 13° K^tt^ decays where 
the K~ and the tt"*" are misidentified as tt^ and , re- 
spectively. Combinatorial background events are those 
not described by the above components; they do not ex- 
hibit any peaking structure in tuxtt or Am. 

The functional forms of the probability density func- 
tions (PDFs) for the signal and background components 
are chosen based on studies of Monte Carlo (MC) sam- 
ples. However, all parameters are determined from two- 
dimensional likelihood fits to data over the full 1.81 < 
rriK^ < 1.92GeV/c2 and 0.14 < Am < O.lGGeV/c^ re- 
gion. 

We fit the RS and WS data samples simultaneously 
with shape parameters describing the signal and random 
7r+ components shared between the two data samples. 
We find 1, 141, 500±1, 200 RS signal events and 4, 030±90 
WS signal events. The dominant background component 
is the random 7r+ background. Projections of the WS 
data and fit are shown in Fig. [1] 
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FIG. 1: a 

0.1445 < Am < 0.1465 GeV/c^ 
dates with 1.843 < ruKn < 1.883 GeV/c^ The fitted PDFs 
are overlaid. The shaded regions represent the different back- 
ground components. 



The measured proper-time distribution for the RS sig- 
nal is described by an exponential function convolved 
with a resolution function whose parameters are deter- 
mined by the fit to the data. The resolution function 
is the sum of three Gaussians with widths proportional 



to the estimated event- by- event proper-time uncertainty 
6t. The random background is described by the same 
proper-time distribution as signal events, since the slow 
pion has little weight in the vertex fit. The proper-time 
distribution of the combinatorial background is described 
by a sum of two Gaussians, one of which has a power-law 
tail to account for a small long-lived component. The 
combinatorial background and real D° decays have dif- 
ferent St distributions, as determined from data using a 
background-subtraction technique [il'l based on the fit to 
rriK-K and Am. 

The fit to the RS proper-time distribution is performed 
over all events in the full m^^ ^ and Am region. The PDFs 
for signal and background in m^^ and Am are used in 
the proper-time fit with all parameters fixed to their pre- 
viously determined values. The fitted lifetime is found 
to be consistent with the world-average lifetime [2^ . 

The measured proper-time distribution for the WS sig- 
nal is modeled by Eq. ^ convolved with the resolution 
function determined in the RS proper-time fit. The ran- 
dom TTj^ and misreconstructed D° backgrounds are de- 
scribed by the RS signal proper-time distribution since 
they are real D'^ decays. The proper-time distribution 
for WS data is shown in Fig. [2l The fit results with and 
without mixing are shown as the overlaid curves. 
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Random Tig 
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FIG. 2: a) The proper-time distribution of combined D'^ 
and Z)° WS candidates in the signal region 1.843 < mx-n- < 
1.883 GeV/c^ and 0.1445 < Am < 0.1465 GeV/c^ The result 
of the fit allowing (not allowing) mixing but not CP violation 
is overlaid as a solid (dashed) curve. Background compo- 
nents are shown as shaded regions, b) The points represent 
the difference between the data and the no-mixing fit. The 
solid curve shows the difference between fits with and without 
mixing. 

The fit with mixing provides a substantially better de- 
scription of the data than the fit with no mixing. The 
significance of the mixing signal is evaluated based on 



the change in negative log Hkehhood with respect to the 
minimum. Figure [3] shows confidence-level (CL) contours 
calculated from the change in log likelihood (— 2Aln£) 
in two dimensions (x'^ and y') with systematic uncer- 
tainties included. The likelihood maximum is at the un- 
physical value of x'^ = -2.2 x 10"^ and y' = 9.7 x lO'^. 
The value of — 2Aln£ at the most likely point in the 
physically allowed region (x'^ = and y' = 6.4 x IQ-^) 
is 0.7 units. The value of — 2Aln£ for no-mixing is 
23.9 units. Including the systematic uncertainties, this 
corresponds to a significance equivalent to 3.9 standard 
deviations (1 — CL = 1 x 10"'') and thus constitutes 
evidence for mixing. The fitted values of the mixing pa- 
rameters and i?D are listed in Table HI The correlation 
coefficient between the a;'^ and y' parameters is —0.94. 




-0.5 



0.0 

x'2/10'^ 



0.5 



FIG. 3: The central value (point) and confidence-level (CL) 
contours for 1 - CL = 0.317 (la), 4.55 x 10"^ {2a), 2.70 x 
10"^ {3a), 6.33 x 10"^ (4cr) and 5.73 x 10"^ {5a), calculated 
from the change in the value of —2ln£, compared with its 
value at the minimum. Systematic uncertainties are included. 
The no- mixing point is shown as a plus sign (-I-). 



Allowing for the possibility of CP violation, we cal- 

^ val- 



and 



R^)/{R^ + R^) listed in TablelH from the fitted 



culate the values of i?D = 

ues. The best fit in each case is more than three standard 
deviations away from the no-mixing hypothesis. All cross 
checks indicate that the high level of agreement between 
the separate D° and 13° fits is a coincidence. 

As a cross-check of the mixing signal, we perform inde- 
pendent {rriKTn Am} fits with no shared parameters for 
intervals in proper time selected to have approximately 
equal numbers of RS candidates. The fitted WS branch- 
ing fractions are shown in Fig. |4] and are seen to increase 
with time. The slope is consistent with the measured 
mixing parameters and inconsistent with the no-mixing 
hypothesis. 

We have validated the fitting procedure on simulated 
data samples using both MC samples with the full de- 
tector simulation and large parameterized MC samples. 
In all cases we have found the fit to be unbiased. As a 
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FIG. 4: The WS branching fractions from independent 
{rriKTv, Am} fits to slices in measured proper time (points). 
The dashed line shows the expected wrong-sign rate as de- 
termined from the mixing fit shown in Fig. [3 The with 
respect to expectation from the mixing fit is 1.5; for the no- 
mixing hypothesis (a constant WS rate), the is 24.0. 



TABLE L Results from the different fits. The first uncertainty 
listed is statistical and the second systematic. 



Fit type 



Parameter Fit Results (/lO 



No CP viol, or mixing 


Ru 


3.53 


± 


0.08 


± 


0.04 


No CP 


Rd 


3.03 


± 


0.16 


± 


0.10 


violation 




-0.22 


± 


0.30 


± 


0.21 




y' 


9.7 


± 


4.4 


± 


3.1 




Rd 


3.03 


± 


0.16 


± 


0.10 


CP 


Ao 


-21 


± 


52 


± 


15 


violation 




-0.24 


± 


0.43 


± 


0.30 


allowed 




9.8 


± 


6.4 


± 


4.5 






-0.20 


± 


0.41 


± 


0.29 




y'- 


9.6 


± 


6.1 


± 


4.3 



further cross-check, we have performed a fit to the RS 
data proper-time distribution allowing for mixing in the 
signal component; the fitted values of the mixing param- 
eters are consistent with no mixing. The correlations 
among parameters determined at different stages of the 
fit are low. In addition we have found the staged fit- 
ting approach to give the same solution and confidence 
regions as a simultaneous fit in which all parameters are 
allowed to vary. 

In evaluating systematic uncertainties in i?D and the 
mixing parameters we have considered variations in the 
fit model and in the selection criteria. We have also con- 
sidered alternative forms of the mx-m Am, proper time, 
and 5t PDFs. We varied the t and 5t requirements. In 
addition, we considered variations that keep or reject all 
candidates sharing tracks with other candidates. 

For each source of systematic error, we compute 
the significance sf — 2 ln£(a;'^, y') — lnC{x''^,y^) /2.3, 

where {x''^, y') are the parameters obtained from the 
standard fit, {x' ^ , y[) the parameters from the fit includ- 
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ing the i systematic variation, and C the likelihood of 
the standard fit. The factor 2.3 is the 68% confidence 
level for 2 degrees of freedom. To estimate the signifi- 
cance of our results in (.t'^, y'), wc reduce — 2Aln£ by a 
factor of 1 + i7s? = 1.3 to account for systematic errors. 
The largest contribution to this factor, 0.06, is due to 
uncertainty in modeling the long decay time component 
from other D decays in the signal region. The second 
largest component, 0.05, is due to the presence of a non- 
zero mean in the proper time signal resolution PDF. The 
mean value is determined in the RS proper time fit to be 
3.6 fs and is due to small misalignments in the detector. 
The error of 15 x 10~^ on Au is primarily due to uncer- 
tainties in modeling the differences between and K~ 
absorption in the detector. 

We have presented evidence for D'^-D'^ mixing. Our 
result is inconsistent with the no-mixing hypothesis at a 
significance of 3.9 standard deviations. We measure y' = 
[9.7 ±4.4 (stat.) ±3.1 (syst.)] x 10"■^ while x'^ is consis- 
tent with zero. We find no evidence for CP violation and 
measure Ru to be [0.303±0.016 (stat.) ±0.010 (syst.)]%. 
The result is consistent with SM estimates for mixing. 
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